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Abstract 

We calculate the excitation of low frequency gravity waves by turbulent convection in 
the sun and the effect of the angular momentum carried by these waves on the rotation 
profile of the sun's radiative interior. We find that the gravity waves generated by con- 
vection in the sun provide a very efficient means of coupling the rotation in the radiative 
interior to that of the convection zone. In a differentially rotating star, waves of differ- 
ent azimuthal number have their frequencies in the local rest frame of the star Doppler 
shifted by different amounts. This leads to a difference in their local dissipation rate and 
hence a redistribution of angular momentum in the star. We find that the time scale for 
establishing uniform rotation throughout much of the radiative interior of the sun is ~ 10 7 
years, which provides a possible explanation for the helioseismic observations that the solar 
interior is rotating as a solid body. 
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1. Introduction 



Our understanding of the internal structure of the sun has improved significantly 
since the advent of helioseismic mapping. The new dedicated observing programs such 
as GONG and SOHO are providing a map of the sun's internal structure with increasing 
resolution and accuracy (Christensen-Dalsgaard et al. 1996). We already have an excellent 
measurement of the thickness of the convection zone (Christensen-Dalsgaard et al. 1996), 
and the rotation profile in the outer half of the sun (Thompson et al. 1996). These results 
suggest that the radiative interior of the sun is rotating as a solid body at a rate equal to 
the rotation rate of the base of the convection zone. It has been suggested that magnetic 
torques in the radiative interior of the sun might be responsible for its solid body rotation. 
In this paper, we show that gravity waves generated by turbulent stresses in the convection 
zone carry significant angular momentum and are very efficient in establishing solid body 
rotation throughout much of the radiative interior. 

The calculation of the mechanical energy luminosity in gravity waves due to turbulent 
stresses in the convection zone is discussed in the next section. The angular momentum 
luminosity in gravity waves and its effect on the rotation profile of the solar interior is 
discussed in section 3. 

2. Energy luminosity in low frequency gravity waves due to turbulent excitation 

We follow the work of Goldreich, Murray, and Kumar (1994) in order to calculate 
the energy luminosity in gravity waves due to excitation by the fluctuating Reynold's 
stress. For simplicity, we assume that the wave excitation is independent of the azimuthal 
wavenumber. The mechanical energy luminosity in gravity waves per unit frequency, i.e., 
the energy flowing through a sphere per unit time and frequency, just below the convection 
zone is given by 



where £ r and £/,, are the radial and horizontal displacement wavefunctions which are nor- 
malized to unit energy flux just outside the convection zone, t is the spherical harmonic 
degree of the wave, v is the convective velocity, L is the radial size of an energy bearing 
turbulent eddy, T L/v is the characteristic convective time, and ( is the ratio of the 
horizontal and vertical length scales of the turbulent eddies. 

For turbulence in the solar convection zone, L ~ aH and ( ~ 1, where a ~ 1.65 
is the mixing length parameter and H is the pressure scale height. Gravity waves of 
long horizontal wavelength, i.e. low £, and frequency to are most efficiently excited by 
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energy bearing turbulent eddies which have turnover times comparable to cu _1 . The energy 
luminosity per unit frequency, as seen by an observer corotating with the base of the 
convection zone, at frequencies greater than the characteristic convective frequency at the 
base of the convection zone (N c fx 0.15/uHz), is readily estimated from equation (1): 



W>~ ' , " {^) [fj ^{^) «* (2) 

where L Q « 4x 10 33 erg s -1 is the solar luminosity and the subscript uj indicates a quantity 
evaluated at r such that uT(r) ~ ir/2. The latter equality in equation (2) is only applicable 
to dipole waves and follows from using r u fx r c fx 0.72R & and H u fx H c (uj/N c )~ 3 /( n+3 \ 
where H c fx 0.1r c is the pressure scale height at the base of the convection zone and n fx 1.5 
is the effective polytropic index of the solar convection zone. This is in good agreement 
with the numerical evaluation of equation (1) which yields a power law of index —4.5. For 

(c) 

u) fx N c , our estimate of L E is, to within a factor of a few, equivalent to that of Press 
(1981), who considered the effect of low frequency gravity waves on elemental mixing in 
the radiative interior of the sun. For uj > N c , however, Press underestimated the energy 
input into gravity waves by only considering the excitation due to inertial range eddies at 
r fx r c ; the Kolmogorov scaling gives oc a; -8 5 for inertial range eddies. 1 

3. Angular momentum transport by gravity waves 

The angular momentum luminosity per unit frequency (angular momentum flowing 
across a sphere of radius r per unit time and frequency), L' a (uj, £, m, r), associated with 
a gravity wave of frequency u and azimuthal order m can be shown to be equal to (e.g., 
Goldreich & Nicholson 1989) 

L^,e,m,r)= - mL ^- e - m - r \ (3) 

where Le is the mechanical energy luminosity per unit frequency associated with the 
wave. It should be pointed out that Le and u are frame dependent quantities, but the 
ratio Le/oj, the luminosity of wave-action, is frame independent. Thus the value of L' a 
does not vary as the wave propagates through a non-dissipative shearing medium. We find 



1 Waves with u> > iV c also get excited by turbulence in the convective overshoot layer. 
If the thickness of the overshoot layer is ~ r]H c then, from equation (1), we estimate that 
the luminosity generated at u> ~ by turbulence in the convective overshoot layer is 

~ ?7~ 4 / 3 times larger than that generated at the same frequency in the convection zone. 
We note that the singularity as r\ — > does not apply to stars since uj for gravity waves 
has an upper bound of order the maximum Brunt- Vaisala frequency in the star. 



3 



it convenient to evaluate the frame dependent quantities, L# & uj, as seen by an observer 
corotating with the base of the convection zone (where the waves are generated). The 
net angular momentum luminosity in dipole gravity waves of a non-zero m at the base of 
the convection zone is obtained by combining equations (2) and (3) and is ~ 10 35 dyne 
cm. This luminosity, deposited in the sun over ps 10 7 years, is equal to the total angular 
momentum in the solar interior (provided that waves of opposite m values do not deposit 
their angular momentum at the same location, which we show below is the case). Thus 
gravity waves should have a significant effect on the interior rotation profile of the sun. 

Differential rotation along latitudes in radiative stars is subject to instabilities which 
restore uniform rotation (e.g. Balbus & Hawley 1996) and so we consider the angular 
rotation frequency in the radiative interior of the sun, O, to be a function of r alone. The 
rotation speed at the base of the convection zone is denoted by fi c . The wave frequency in 
the rest frame of the fluid at radius r, u;* (r), is given below in the terms of its frequency 
u at the base of the convection zone 

U m) (r) = u + m[Q c - Q(r)]. (4) 



If we ignore wave dissipation and assume that there are no corotation resonances in the 
system, i.e., u>i m \r) is non-zero everywhere, then the net angular momentum luminosity 
carried by gravity waves is zero everywhere in the radiative interior and therefore waves 
do not modify the rotation profile of the sun. 

However, the radiative dissipation rate of gravity waves depends on their frequency 
in the local rest frame of the fluid, which is different for waves of ±m in a differentially 
rotating star. This leads to a non-zero net angular momentum luminosity whenever the 
rotation profile of the star deviates from solid body rotation. 

The local radiative dissipation rate for gravity waves can be expressed as 

(u i r ) » ^ ( dlnT \ (9lnT\ _ 2F r k 2 H T 
^ ' ' pc 2 s dlnT/dr \dlnp J \dlnp ) s hp 

where k r « Ny/£(£ + l)/[rivi m) (r)] is the wave's radial wave number in the local rest frame 
of the fluid, F r is the radiative flux a distance r from the center, c s is the sound speed, p 
is the pressure and Ht is the temperature scale height. 

From equation (5) we see that low frequency gravity waves with u ~ N c are strongly 
damped near the base of the solar convection zone. For the solar model we are using (due 
to Christensen-Dalsgaard) the damping length of dipole waves with uj ps N c is only ~ 10 
wavelengths. Waves with uj ~> N c carry less energy, but their damping length is longer 
and so they propagate deeper into the radiative interior. The energy luminosity per unit 
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frequency in gravity waves of frequency u, degree £, and azimuthal order m, a distance r 
from the center, is given by: 

Le(uj, £, m, r) = L^\u>, £) exp[— r(u, £, m, r)], (6) 

where 

f r ° <y(uji m) (r') £ r') 
r(u,e,rn,r)= / dr' ^J^/^ (7) 



is the attenuation depth of the wave and v g ~ u)i m \r)/k r is the group speed of gravity 
waves in the radial direction. We note that ^/v g oc uj~ 4 [£(£ + l)] 3 / 2 . The net energy 
luminosity is obtained by integrating Le(uj, £, m, r) over all frequencies and adding the 
contribution from waves of different £ and m. For the case of solid body rotation, the 
net energy luminosity in the sun is shown in Figure 1, which also shows the frequency for 
which r ~ 1 for dipole waves. Since wave excitation falls off with increasing frequency 
and £, and the attenuation depth decreases as ~ £ 3 o;~ 4 , most of the contribution to the 
luminosity comes from dipole gravity waves with frequencies such that r ~ 1. 

The net angular momentum luminosity at r, L a , is obtained by adding waves of 
different m & £, and integrating over wave frequency 



f°° \m\L (yC \u £) ( \ 
L *( r )=J2l duJ u 1 exp[-r(o;,l,-lm|,r)] - exp[-T(u,£,\m\,r)] 1. (8) 

£,\m\ ® \ J 



In the limit of small differential rotation, |AO(r)| = |fi(r) — O c | <C uj, this equation 
simplifies to 



L a(r) ~ V / du 



2\m\L e(oj , £, 0, r) 



sinh(|m|5r), 



(9) 



where 



5T = i r dr ,:(^y)An { r>) 

J r Vg(uj,£,r') lo 



Thus the rate at which angular momentum is deposited in the fluid is 



(10) 



dL a {r) 
dr 



£,\m\ 



\m\i(u,l,r)L E {w,l,Q,r) 
uv g {u,£,r) 



4\m\ ( j cosh(|m|5r) — sinh(|m|5r) 



(11) 



Finally, the rate of change of the rotation frequency is given by the following equation 
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d_ 

~dt 



[Afi(r)] 



3 dL a (r) 
8irpr 4 dr 



(12) 



Equations (11) and (12) together determine the evolution of the rotation profile of the 
sun. The resulting equation is a nonlinear differentio-integral equation that is not easy to 
solve. We have, however, numerically solved the linearized version of this equation, which 
is valid in the limit that St <C 1; the result for the evolution of the angular rotation rate 
of the sun is shown in Figure (2). Since r ~ 1 for the waves with most energy and since 
St < t(AQ/u), St < 1 is in fact a good approximation in the solar interior. 

We now focus on providing simple estimates of the timescale over which the rotation 
rate changes and a physical discussion of the angular momentum redistribution in the sun. 
For small differential rotation (|AO|/u; <C 1) L a (r) increases linearly with AO and so there 
is a characteristic time, t(r), for the change of the angular momentum at radius r that is 
independent of the rotation rate 



Figure (3) shows t(r) for the sun. Note that the characteristic time in the radiative interior 



influenced by the angular momentum deposited by gravity waves. To gain some physical 
understanding of how the rotation profile evolves, we go back to equation (9). 

Note that Le is a decreasing function of (r c — r) (waves are rapidly attenuated with 
distance; see Fig. 1) while sinh(|m|5r) is an increasing function. Thus the angular mo- 
mentum luminosity, which is proportional to the product of these functions (see eq. [9]), 
peaks at some depth that depends on the rotation profile, and so the sign of the angular 
momentum deposited changes at this depth. At small depths (before the peak) the angular 
momentum deposited has a sign such as to reduce |AO|. The angular momentum removed 
from this region gets deposited deeper in the star over a distance which is of order a wave 
damping length (see Fig. 2). 

Initially gravity waves force a shell of thickness equal to a dissipation length of waves 
with uj rs N c , lying just below the convection zone, to corotate with the convection zone; 
this occurs on a time scale of only ~ 10 3 years (fig. 3). The angular momentum removed 
from this region is deposited over the next dissipation length. This process continues and 
the thickness of the shell corotating with the convection zone grows with time. Since 
the dissipation length increases with depth, the angular momentum removed from the 
corotating shell is redistributed over increasingly larger distances. These effects are all seen 
in Figure (2). On a timescale 10 7 years the size of the differentially rotating core of the 
sun has shrunk to less than ps 0.5-Rq. Will this process lead to the entire radiative interior 




(13) 



is less than m 10 



years and so the rotation rate in this region is expected to be strongly 
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rotating uniformly? For simplicity, we have assumed that waves of different m are excited 
to the same amplitude. This means that the net angular momentum luminosity at the base 
of the convection zone is zero, and thus any excess/deficit of angular momentum in the 
radiative interior (over corotation with the base of the convection zone) gets concentrated 
in a region near the center of the sun the size of which decreases with time. Moreover, 
we have not considered the effect of waves which return to the convection zone from the 
radiative interior; these waves will be reabsorbed by the convection zone, thus removing 
the excess/deficit of angular momentum from the core of the sun. Thus, we believe that 
allowing for a net exchange of angular momentum between the convection zone and the 
radiative interior will lead to uniform rotation of the entire radiative interior on the time 
scale of ~ 10 7 years estimated above. 

4. Discussion 

We have calculated the excitation of low frequency gravity waves by turbulent convec- 
tion in the sun and the angular momentum luminosity they carry. The energy luminosity 
in dipole gravity waves with frequencies of order N c ~ 0.15/uHz (the characteristic con- 
vective frequency at the bottom of the convection zone) is ps 10 29 erg s _1 near the base 
of the convection zone. Higher frequency dipole and quadrupole gravity waves are excited 
most efficiently by energy bearing turbulent eddies deeper in the convection zone, and the 
resultant power spectrum is cu -13 / 3 . Waves of higher £ are more evanescent in the convec- 
tion zone and are thus less efficiently excited. Low frequency gravity waves suffer strong 
radiative damping in the solar interior; dipole waves with frequencies of ~ 0.15/uHz travel 
only ~ 10 wavelengths (~ 0.02Rq) into the radiative interior, whereas I = 1(10) waves of 
frequency 1.0/uHz have dissipation lengths ~ 1.0(0.1) Rq. 

The angular momentum luminosity associated with a wave, i.e., the amount of angular 
momentum crossing a sphere per unit time, is m times the ratio of the energy luminosity 
and the wave frequency (where m is the azimuthal number of the wave). The local damping 
rates for ±m waves are in general different in a differentially rotating medium, which 
results in a local deposition of angular momentum. We find that there is enough angular 
momentum in gravity waves generated by convection that they can force much of the 
radiative interior of the sun into corotation with the base of the convection zone in ps 
10 7 years (see Fig. 2 & 3); Zahn et al. (1996) have independently arrived at a similar 
conclusion (see also, Schatzman 1996). 

Our calculations were carried out under the assumption that waves of different az- 
imuthal number are excited to the same amplitude, as seen in the local rest frame of the 
sun just below the convection zone. Furthermore, we have neglected the effect of the cori- 
olis force on the wave function as well as the angular momentum removed from the core of 
the sun by gravity waves which are reabsorbed by the convection zone. As a consequence 
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waves merely redistribute angular momentum in the radiative interior of the sun in such 
a way that a significant fraction of the radiative interior is forced into corotation with 
the bottom of the convection zone, where the waves are excited. We feel that when these 
simplifying assumptions are dropped, and a net exchange of angular momentum between 
the convection zone and the radiative interior is allowed, the mechanism we have described 
here will lead to uniform rotation of the entire radiative interior of the sun. 
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Figure Captions 



FIG. 1. — (a) The net energy luminosity in gravity waves (L E ) as a function of radius in the 
sun. Wave generation in the convective overshoot layer is not included, (b) The frequency 
of dipole waves, as a function of radius in the sun, for which the energy luminosity in the 
wave has decreased by a factor of e from its value at the base of the convection zone (i.e., 
the wave's attenuation depth, r, is ~ 1). These waves carry most of the energy luminosity 
at the given radius. We note that uj t= ± oc [£(£ + 

FIG. 2. — (a) The time evolution of the differential rotation profile of the sun, AO(r), which 
is obtained by solving eqs. (11) & (12). The curve marked with t = is an arbitrarily 
chosen initial rotation profile and is not intended to reflect the sun's actual primordial 
rotation profile. The energy luminosity in gravity waves, needed for this calculation, is 
given in Figure (1). Panel (b) shows the net angular momentum luminosity in gravity 
waves, L a (r), at the three different epochs considered in panel (a). 

FIG. 3. — t(r) is the characteristic time over which the solar rotation rate at radius r 
changes due to the angular momentum deposited by gravity waves. t(r) is much less than 
the age of the sun everywhere in the radiative interior and so the angular momentum 
carried by gravity waves can strongly modify the solar rotation profile. 
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